model {

for(i in 1:N){
mu_M[i] <- alpha*x[i]
M[i] ~ dnorm(mu_M[i],prec1)

logit(mu_y[i]) <- c*x[i]+beta*M[i]
y[i] ~ dbern(mu_y[i])
}

alpha ~ dnorm(0, 1.0E-3)
beta  ~ dnorm(0, 1.0E-3)
IE<-alpha*beta
c ~ dnorm(0, 0.01) 
DE<-c
var1 ~ dgamma(1,0.1)
prec1 <-1/var1
}